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Multiband superconductors support the excitation of vortices with fractional quantum flux, known as frac- 
tional vortices. In the ground state, the fractional vortices in different bands bond together to form a composite 
vortex with the standard flux quantum cDq; thus, it is difficult to stabilize the fractional vortices. Here we show 
that fractional vortex lattices can be stabilized in multiband superconductors with a periodic pinning array at the 
half matching field of the composite system, when full matching of the fractionalized vortices occurs. In the 
presence of a high current, the fractional vortices in different bands decouple and move at different velocities. 
When the current is turned off suddenly, the fractional vortices in different bands may be trapped at different 
pinning sites. This system also exhibits rich dynamic behavior, and for fractional vortices in two-band supercon- 
ductors we find a phases where the vortices one band are pinned and the vortices in the other band are moving. 
These different phases can be observed in transport measurements, and produce a new type of self-induced 
Shapiro steps that can arise under the application of only a dc drive. 
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PACS numbers: 74.25.Ha, 74.25.Uv, 74.70.Ad 

Multiband superconductivity continues to attract enormous 
attention as more and more multiband superconductors are 
discovered, such as MgB2' and iron-based superconductor J^. 
Multiband superconductors have unique features that cannot 
be found in single band superconductor^"'^. Examples in- 
clude fractional vortices^ '■*, short-range repulsion and long- 
range attraction between vortice^^^ oscillations of supercon- 
ducting phase difference between different bands (the Leggett 
mode)"* '^ and the time-reversal symmetry broken phase'^, just 
to name a few. 

In multiband superconductors, electrons in each band form 
a superfluid condensate characterized by a complex gap func- 
tion exp(/0^). It is possible that the phase 0^ of the differ- 
ent condensates can wind by 21 jn with an integer Ij. If the 
phase of the y-th band winds by 2lj7T, it caiTies fractional flux 

~ huT- [Ijf/ ^^/'"/jj *^o, where mj is the electron's effec- 
tive mass and - he lie is the flux quantum^. When the 
winding numbers in different condensates are the same, Ij - I, 
a composite vortex appears carrying ZOq fux. The composite 
vortex can be viewed as a bound state of fractional vortices in 
different bands with the vortex cores located at the same posi- 
tion. For multiband superconductors under a strong magnetic 
field, the ground state is a lattice of composite vortices with 
Ij - 1 . Recently there has been growing interest in the stabi- 
lization of fractional vortices in multiband superconductors. 

It was shown that the energy of the fractional vortices di- 
verges logarithmically or linearly in bulk superconductors^, 
which renders the fractional voitices unstable in bulk. This 
logarithmic divergence is cut off in a finite sample or in a 
mesoscopic superconductor. Thus the fractional vortices can 
be stabilized in mesoscopic superconductors, as demonstrated 
in Refs. ri"8l420l bv numerical solution of the Ginzburg-Landau 
equation. The fractional vortices can also be stabilized near 
the suiface of a multiband superconductoi'^il It was demon- 
strated that under a driving cuiTent, the composite vortex lat- 
tice dissociates if the disparity of the fractional flux and co- 
herence length in the different bands is large^^, resulting in 
flux flow of the decoupled fractional vortex lattices. The com- 



posite vortex lattice can also be decoupled by thermal fluctu- 
ations. As temperature increases, the fractional vortex lattice 
in the band with a weaker phase rigidity first melts while the 
fractional vortex lattice order in the band with a stronger phase 
rigidity remainJ^ 

In this work, we show that fractional vortices can be stabi- 
lized in bulk multiband superconductors with a periodic array 
of pinning sites. For square periodic pinning arrays, peaks in 
the critical current occur when the number of vortices is an 
integer multiple of the number of pinning siteP^EHl as well as 
at fractional fields such as 1 /4, 1 /2 and 3/4, where crystalline 
vortex states form that can match with the symmetry of the 
pinning arrajEEl We consider a square pinning array with 
matching field under an applied field creating a composite 
vortex lattice with B - B^/2, so that when the vortices are 
fractionalized, the total number of vortices will equal the first 
matching field. 

To stabilize a fractional vortex lattice we apply a current 
large enough to cause the fractional vortex lattices in different 
bands to move at different velocities. When we suddenly turn 
off the current, the fractional voitices in different bands can 
be trapped by different pinning sites. Although the resulting 
fractional vortex lattice is metastable, the lifetime is long if the 
pinning is strong and the final vortex number equals the first 
matching field. In addition, we show that the dynamic phase 
diagram for voitices in multiband superconductors exhibits a 
wide range of dynamical behaviors, including a phase where 
the vortices in one band are pinned while the voitices in the 
other band are moving. We found that this system can exhibit 
novel self-induced Shapiro steps in the flux flow of fractional 
vortex lattices with different velocities, even though only a dc 
current is applied. 

We consider a two-band superconductor with interband 
Josephson coupling. To demonstrate our idea we adopt the 
phenomenological London theory, which is valid when A is 
much larger than the coherence length in different bands, 
such that the amplitude of the order parameter is approxi- 
mately constant in space. The London free energy functional 
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where X,, 



^(}n^c'^~)l (47rn^e2) is the London penetration 
depth for each condensate with superfluid density n^. A is the 
vector potential and y is the interband Josephson coupUng. 
The effective penetration depth for the two-band system is 
X'^ - 2^=1 /i^^. For MgB2, y > 0; thus, in the ground state, 
6*1 - 02- For certain iron-base d sup erconductors, y < 0; thus, 
in the ground state 0\ = 02 + the proposed liquid hy- 

drogen superconductorjIS^ISlthe interband Josephson coupling 
is absent"*^, y = 0. 

In the vortex lattice phase, the profile of the superconduct- 
ing phase is given by the following equation obtained by 
minimizing Eq. ([T]i with respect to 0^, 



02)-rsin(6li -02) = 0, 



(2) 



subject to the boundary condition V x (V0^) = In Yj^j S{r - 
Yfij), where = A^<^()/A^ is the fractional quantum flux and 
Tfij is the core position of the fractional vortex in the band 
fi. We assume that the vortex lines are straight. At higher 
magnetic fields the Josephson coupling becomes negligible'^. 
Neglecting the Josephson coupling y = 0, Eq. Q can be 
rewritten afi^ 



+ aHV X Bf] + [V (01 - 02)f . 



(3) 



The energy only involves quadratic terms; thus, Eq. Q de- 
scribes the pairwise interaction between vortices. The first 
term is the electromagnetic coupling between vortices and the 
second term accounts for the coupling between the supercon- 
ducting phases. The electromagnetic coupling contributes to 
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FIG. 2. (color online) (a) Dynamic phase diagram for two-band su- 
perconductors with a pinning array subject to a dc current. Four dif- 
ferent regions can be identified. Here = and 02/'l'i = 5.0. 



the repulsive interaction between fractional vortices in differ- 
ent bands and in the same band. The term proportional to 
(V0^)^ describes the repulsive interaction between vortices in 
the same band, and the cross term proportional to -'V6i'V02 
describes the attraction between vortices in different bands. 
This attraction outweighs the repulsion due to the magnetic 
coupling and the net interaction between vortices in different 
condensates is attractive. The repulsion between the fractional 
vortices in the same band at a separation ^ = Tf^j - r^j 
thus can be written a^3^ 
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and the attraction between two vortices in the different con- 
densates with a distance ri2,,y = rj ,■ - r2j is 
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where KqH) is the modified Bessel function of the first kind. 
Equations Q and (jSjl are valid away from the vortex core. 

We then introduce pinning arrays into the system. The in- 
teraction between the pinning site at Vp and the fractional vor- 
tex in the //-th band at is modeled as 



(r^, - r^) = -r^ exp -(r^ - Vp]^ III 



(6) 



FIG. 1. (color online) (a) Dependence of the velocity of fractional 
vortices on current and (b) IV characteristics for different pinning 
strengths F^. For a weak pinning strength, the two vortex lattices de- 
pin simultaneously at the depinning current and then move together 
with the same velocity. For a large current, they split into two mov- 
ing fractional vortex lattices with distinct velocities. In the case with 
a stronger pinning, they depin at different currents. Here jip = and 
(D2/<l'i = 5.0. 



where characterizes the pinning strength and is the pin- 
ning range. The equation of motion for the fractional vor- 
tices in different bands thus is governed by the following over- 
damped dynamics 



(7) 



where the last term on the right-hand side of Eq. (j7]i is the 
Lorentz force from a current density J. The fractional vortices 
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FIG. 3. (color online) (a) Current-velocity curves and (b) the cor- 
responding differential resistivity dEldJ for different (5)2l'^i- The 
spikes in the velocity are due to the resonance among the two mov- 
ing fractional vortex lattices and the pinning array. The first two 
peaks from the left in (b) are due to the depinning transition of the 
two composite vortex lattices. The other peaks in (b) are due to the 
Shapiro steps. Here np = iiy and Ti = = 2.0. 



with a larger flux experience a stronger Lorentz force. 77^ is 
the viscosity of the fractional vortices. The viscosity depends 
on the size of normal core and is given by 77^ = <t>Q/(27Tc^^J^) 
following the Bardeen-Stephen model^^. In the absence of a 
pinning potential f/^ p - 0, it was shown in Ref. |22jthat the 
two fractional vortex lattices move with the same velocity for 
a low current. At a high current, the two fractional vortex lat- 
tices decouple from each other resulting in a decoupled flux 
flow of fractional vortex lattice. In the presence of an ac cur- 
rent supeiposed on a dc cuiTent, Shapiro steps are induced. 

We solve Eq. (|7]i numerically to obtain the velocity of frac- 
tional voitices at a given J. From the power balance condi- 
tion 771 + 772 = JEa^ with E the electric field and a the 
vortex lattice constant, we obtain the electric field. We in- 
troduce dimensionless units for force: <I>i<l>2/(87r^/l"'); length: 
A; time: 87r^77i/i'^/(<l)i<l)2); and current: c<l>2/{^n-A^), and we 
set 771 - 772 and = 1. We consider two pinning densities 
Hp = Til, and Hp - Itiy, where riy is the density of the compos- 
ite vortices. In the case Hp — we use a square pinning array 
with a lattice constant Up^ - Op^y = 5.0. For Hp - Iriy, we 
use a rectangular pinning array with y = 2flp_v - 5.0. The 
total number of composite vortices is 256. We use a square 
simulation box with size 80 x 80 and apply periodic boundary 
conditions. Equation (j?]) is integrated using the second order 
Runge-Kutta method with a time step A? - 0.02. 

The attraction between two fractional vortex lattices in 
different bands is a periodic function of space with a pe- 
riod equal to the lattice constant. The maximal attraction 
is Fa = <l)i<l)2a/(647r*'/l'^).l^ As shown in Fig. [l] obtained 
numerically, in the case when the maximal pinning force 
Ff^ p - -'^U^^p is much smaller than F^, i.e. ^ <K F^, the 
two fractional vortex lattices depin simultaneously at a current 
Jd = {Fi,p+F2,p)c/<I>o. After depinning, they move with equal 
velocities until the current is large enough to decouple them 
as in the case of clean systems. When they are both pinned 
or move with the same velocity, they bond together to form 



a composite vortex lattice. Once they move at different ve- 
locities, the composite vortex lattice splits into two fractional 
vortex lattices. In the IV curve the voltage is zero in the pinned 
phase if one neglects the creep motion of voitices. A nonzero 
electric field is induced when the fractional vortex lattices 
move at the same velocity. When the composite vortex lattice 
splits, the diff'erential resistivity increases, as shown in Fig. [T| 
(b). In the other limit when ^ » Fa, the two fractional 
vortex lattices depin at different currents Jp d - Fp^pc/f^f^, and 
they move at different velocities once depinned. The dynamic 
phase diagram is constructed in Fig. |2] As increases from 
0, the depinning cuiTent increases. Meanwhile the region of 
flux flow of the composite vortex lattice shrinks and finally 
disappears. Then the two fractional vortex lattices depin sep- 
arately. At large current 7 » 1, the two fractional vortex 
lattices move at different velocities ss J^^plicrj^). All these 
dynamical phase transitions appear as features in the TV curve. 

When both fractional vortex lattices in bands 1 and 2 are 
moving with distinct velocities at a large J, the slower mov- 
ing lattice in band 1 experiences a periodic potential induced 
by the pinning array. It also feels a periodic force due to the 
fast moving lattice in band 2. Resonance between these two 
periodic drivings occurs when the periods satisfy liT\ - l2T2- 
Here Ti = a/vi is the peiiod of the ac force due to the pinning 
array acting on the vortex in band 1, and T2 - a/{v2 - vi) is 
the period of the ac force due to the fast moving lattice in band 
2. We have assumed that the lattice constants of the pinning 
array and the vortex lattice are equal. In terms of Vp the res- 
onance condition becomes (I2 + li)vi = /2V2- The resonance 
manifests as self-induced Shapiro steps in the TV curve. In the 
numerical simulations, we observe the resonance as shown in 
Fig. [3] In velocity or electric field curves, we observe spikes 
as a consequence of the resonance whenever Vp satisfies the 
resonance condition. The resonance can be seen more clearly 
in the differential resistivity dE/dJ shown in Fig. |3](b). These 
self-induced Shapiro steps are distinct from the typically ob- 
served Shapiro steps'*" in that they occur without the appli- 
cation of an external ac drive. The presence of self-induced 
Shapiro steps is thus an experimentally observable signature 
of the decoupling transition between two vortex lattices. 

So far we have focused on the commensurate pinning case 
Hp - iiy, where the pinned bound vortex lattice forms a square 
arrangement. To optimally stabilize fractional vortices, we 
consider the case Up - Itiy, where in the ground state the 
composite vortices are located in a checkerboard pattern with 
every other pinning site filled, as shown in Fig. 4(a). We ap- 
ply a cuiTent large enough to cause the fractional vortex lat- 
tices in bands 1 and 2 to move at different velocities. We 
then suddenly turn off the current, and find that most of the 
fractional vortex in bands 1 and 2 become trapped at differ- 
ent pinning sites, forming a state in which every pinning site 
is occupied by a vortex. Thus, another square lattice forms 
that now has one-to-one matching with the pinning array. In 
this way, one may realize static fractional vortices that could 
be measured experimentally using various imaging techniques 
such as a SQUID. The resulting fractional vortices are in a 
metastable state; however, the lifetime can be long for strong 
pinning, and is further lengthened due to the commensuration 
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FIG. 4. (color online) (a) Schematic view of the periodic current 
drive used in the simulations. Snapshot of the vortex configuration 
(b) in the ground state and (c) after a current quench. Open circles 
are pinning sites; blue and red circles represent the fractional vortices 
in difi"erent bands. Some fractional vortices are trapped at different 
pinning sites after the current quench. Here ;ip = 2«,,, Fi = r2 = 2.0, 
= = 0.5, Ji = 3.2, Oz/Oi = 5.0 and r = 100. 



between the vortex configuration and the pinning sites, which 
reduces the effects of fluctuations from vortex-vortex interac- 
tions. In simulations, we apply a periodic square wave current 
with 7(f) - Ji forO < t < T and J{t) = for t < f < 2r. When 
7 = 0, we take a snapshot of the system, as shown in Fig.4(b), 
which indicates that most of the fractional vortices are trapped 
at different pinning sites. We have also performed simulations 



with other fillings such as rip/n,, - 3 and rip/ny - 4, or using 
a periodic trapezoid wave form for the driving current J. All 
the results are qualitatively the same as those in Fig. [4] 

Let us discuss the possible application of the present theory 
to the stabilization of a half-quantum vortex in the p-wave su- 
perconductors, which have the potential for topological quan- 
tum computing because they host the Majorana fermion. For a 
small ratio of the spin superfluid density to superfluid density 
Psp/Ps < 1, it was shown that the half-quantum vortex can be 
stabilized in a mesoscopic p-wave superconductoi^. Later 
this half-quantum vortex was observed expeiimentall}'^ in 
Sr2Ru04, which is a strong p-wave superconductor candidate. 
It was then proposed that half-quantum vortices may be stabi- 
lized in the ground state of perforated superconducting films, 
which can be observed by magnetoresistance measurement'*'^ 
For a large Psp/ps > 1 where the two-half quantum voitices 
form a composite vortex by overlapping their normal cores, it 
is difficult to split them by the application of an external cur- 
rent as proposed in this work, because the two-half quantum 
vortices have the same flux and coherence length, and thus 
they respond to the external current in the same way. 

The creation of fractional voitices replies on the splitting 
of the composite vortex, which is possible only for a weak 
interband coupling. Weakly coupled multiband superconduc- 
tors were found to be realized in V3Si'*'* and FeSei_J^, and 
also possibly in MgB2.'^The external field Ha where the in- 
terband Josephson coupling becomes negligible is given by 
the condition V3W^4 Aj = ^/^i^^J(16n^^i^y\). For 
MgB2, it requires Hg > 4 T at temperature T - OK for typical 
parameters. The cuiTent where the composite vortex lattice 
dissociates is Jj ~ 5 x 10'' A/m^, which is much smaller than 
the depaiiing cuiTent. For VsSi and FeSei_x, the required field 
and the dissociation current are smaller To stabilize 
fractional vortices in the pinning sites, the density of pinning 
sites should be comparable to that of the composite vortices. 
For MgB2, the period of pinning array is about 20 nm. 

To summarize, we propose a method to stabilize fractional 
vortex lattices in multi-band superconductors with periodic 
pinning arrays. Under an applied drive, we identify sev- 
eral distinct dynamic phases that occur, including a phase 
where vortices in both bands are pinned, a phase in which 
one species is pinned and the other is moving, and a higher 
drive phase where both species are in motion but at different 
velocities. The different depinning transitions can be observed 
in transport measures such as dV/dl curves. Another feature 
of the moving state is that this system exhibits self-induced 
Shapiro steps in the transport curves that can arise without the 
application of an external ac drive. A fractional vortex lattice 
can be stabilized when the composite system is initially at the 
half-matching field of the periodic pinning aiTay, so that only 
every other pinning site captures a composite vortex. Once the 
system is in the moving state, the drive is quenched to zero and 
a fractional vortex lattice appears in which every pinning site 
captures a fractional vortex. 
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